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ABSTRACT 



We present new Chandra observations of a high redshift galaxy cluster 

discovered in the Red-Sequence Cluster Survey (RCS): RCS043938-2904.7. X- 
ray luminosity measurements and mass estimates are consistent with Lx-Tx and 
Mg-Tx relationships obtained from low-redshift data. Assuming a single cluster, 
X-ray mass estimates are a factor of ~ 10-1 00 below the red-sequence optical 
richness mass estimate. Optical spectroscopy reveals that this cluster comprises 
two components which are close enough to perhaps be physically associated. We 
present simple modeling of this two-component system which then yields an X-ray 
mass and optical richness consistent with expectations from statistical samples 
of lower redshift clusters. An unexpectedly high gas mass fraction is measured 
assuming a single cluster, which independently supports this interpretation. Ad- 
ditional observations will be necessary to confirm the excess gas mass fraction 
and to constrain the mass distribution. 

Subject headings: galaxies: clusters: individual (RCS043938-2904.7) — X-rays: 
galaxies: clusters 
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Introduction 



There is considerable interest in measuring the evolution of the galaxy cluster mass 
function with redshift. In principle, such a measurement can test structure formation and 
constrain the cosmic expansion history. One of the observational challenges in such an 
undertaking is to identify all clusters in a particular volume; the other is to reliably infer 
cluster mass from observable cluster properties, over a range of cluster redshifts. R ecent 
optical surveys, such as the Red-S equence Cluster Surv ey (RCS. iGladders &: Yeell2005l ): and 
the Sloan Digital Sky Survey (e.g.. iKoester et al.ll2007l ) have found large numbers of clusters 
at moderate to high redshift, with 0.2 ^ z < 1 and < z < 0.3, respectively. At the 
high redshift extent, where spectroscopic information can be difficult to obtain on a survey 
scale, various optical properties, such as richness or optical luminosity, can be taken to be 
mass proxies, and have been shown to be strongly correlated with observables understood 
to b e related to cluster raass, such as ga l axy velocity dispersion and X-ray te mperature 



[e.g.. lYee fc Ellingsonll2003l : iLin et al.ll2003l : IPopesso et al.ll2004l : iLopes et al.ll2006r ). In order 
to use optical richness (or any other optical property) as a mass proxy, understanding the 
evolution of the mass-richness relation is critical to extracting cosmological information, since 
any evolution in the relation must be separated from the evolution of the mass function. In 
particular, some care must be taken when attempting to extend the mass/observable relation 
to high redshift as many of these relations are computed using local data. 

Previous work has shown that optically-selected, high-redshift clusters at a fixed optical 
richness are t ypically lower in temperature and underluminous in X-rays when compared t o 
local clusters Jflicks et al.lboosi : Icilbank et aPbood : llubin et allbooi IPonahue et aPboOlj ) . 
Since these clusters generally follow the local luminosity-temperature relationship, the low 
temperatures and luminosities are interpreted as an indication that the cluster galaxies are 
not in virial equilibrium with the X-ray emitting gas. However, it is important to account 
for the significant scatter in the relation between optical richness and X-ray luminosity, as 
an Eddington-like bias will sc atter points of lowe r than expected X-ray luminosity into an 
optically-selected sample (e.g.. lBower et al.lll994j ). Observing these high redshift clusters is 
thus an avenue to study how clusters approach equilibrium. 

RCS043938-2904.7 is part of a larger sample of high redshift, high richness RCS clusters 
selected for X-ray follow up observat ions. X-ray prope rties of the sample as a whole are 
described in a forthcoming paper by iHicks et al.l (120071 ). Here we present Chandra X-ray 
observations and optical spectroscopy of this apparently rich. X-ray underluminous object. 
In §2 we describe the X-ray data reduction and analysis. In §3 we compare mass estimates 
from the X-ray data to the known mass-temperature and mass-richness relations in literature. 
In §4 we describe the optical spectroscopy, in §5 we discuss our results and we summarize 
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in §6. Throughout the paper we take a standard ACDM cosmology with Hq = 70 km s~^ 
Mpc~^, Qrn = 0.3 and = 0.7. At the cluster redshift, z = 0.9558, 1" corresponds to 
a metric distance of 6.64 kpc. Errors are indicated with 90% confidence intervals unless 
otherwise noted. 



2. X-ray observations & data analysis 
2.1. Reduction 

The X-ray data consists of two Chandra ACIS-S3 observations. A log of the observations 
of RCS043938-2904.7 is presented in Table [H Both observations were performed in the 
VFAINT mode, and the CTI correction was applied with a focal plane temperature of - 
120°C. Data reduction was done using CIAO version 3.2 tools and CALDB 3.3.0 starting 
with level 2 event lists from the standard Chandra pipeline. 

The first observation of RCS043938-2904.7 (obsid 3577) had a highly variable back- 
ground and the second observation was kindly provided by the Director of the Chandra 
X-ray Center to compensate. This second observation (obsid 4438) had a low background 
and was analyzed by standard lightcurve methods to select the acceptable background rate 
intervals, keeping 28.7 ks of the total 28.8 ks. More care had to be taken with obsid 3577 
because standard lightcurve methods eliminated more than three-quarters of the data as 
unusable. 

In order to maximize the useful portion of the data, we took the following approach. We 
binned the data into 100-second segments and sorted them from lowest number of counts to 
highest number of counts. We assume that the cluster fiux is constant in time. This means 
that the signal-to-noise ratio is highest in the bin with the least number of counts and lowest 
in the bin with the most counts. Suppose then that there are S signal counts in each bin 
and On total counts in the first n bins. Since On monotonically increases with n, we want 
to find the lowest value of n for which 

Sn S(n + 1) 



V Cn-\ 



1+1 



or equivalently. 



^ > ' (2) 



Cn 



meaning that adding the data in the (n + 1)^^ bin will reduce the total signal-to-noise ratio 
in our filtered data set. This determined the time intervals we included in our final data set 
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independent of the source count rate. In total, 64.6 ks of the total 76.2 ks of data (85%) 
from obsid 3577 were kept. These data were combined with the second observation (obsid 
4438) using the CIAO script merge_all for a total good-data exposure time of 93.3 ks. The 
same script was used to create an exposure map and exposure-corrected image. 

The X-ray data were then filtered to the energy range of 0.3-7 keV and point sources 
were identified and removed using wavdetect. Pulse-height spectra and associated response 
files were created from a 100 pixel (49.2") radius circular region around the X-ray centroid. 
Background spectra were obtained from another 100 pixel radius circular region over 400 
pixels (196.8") away from the cluster center, though still from the same merged image. 
Care was taken to choose the background region in a location where the t wo observation s 



overlapped. Spectral analysis was performed using XSPEC version 12.2 .0 fjArnaud 



and cluster emission was modeled as an optically-thin thermal plasma (IMewe et al. 



1996h 



11985 



19861 ). The redshift was fixed at the value determined from op tical spectroscopy and th e 
heavy element abundance was fixed at 0.3 times Solar abundance (lAnders fc Grevesselll989l ). 
An attempt was made to determine the metal abundances, but the data do not provide a 
significant constraint. The spectral model also included Galac tic photoabsorption. The nH 
FTOOL from NASA's HEASARC Jpickev fc LockmanI ll99ol ) was used to fix the column 
density of HI. Spectra were binned so that a minimum of 20 counts were in each energy bin. 



2.2. Temperature, surface brightness, mass profile & gas mass fractions 

The best fit temperature, fluxes and luminosities (with 90% confidence intervals) are 
listed in Table [2J The bolometric luminosity was estimated by integrating the model spectra 
between 0.05 and 50 keV, and over the 49.2" radius aperture. The corresponding metric 
aperture radius at the cluster redshift is 391 kpc for RCS043938-2904.7. 

We extracted a radial profile of the cluster emission from the exposur e-corrected image. 



Point sources were removed using the CIAO script dmf ilth and a /9-model ( ICavaliere fc Fusco-Femiano 



19761 ) was fit to the data. Data for the fit was taken from a circular aperture with a radius 
of 98.4". This aperture was split into circular annuli with 4.92" thicknesses and fit using the 
Sherpa program included in the CIAO software The /3-model for the surface brightness is: 

3/3+ i 

/(r)=/oll+(-) 1 . (3) 




We used the observed flux within the 49.2" radius spectral aperture to express the normaliza- 
tion in terms of physical units. We also determine the central electron density, n^o, using the 
thermal-plasma emission model discussed above and assuming that the plasma is isothermal 
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throughout the intra-cluster medium (ICM). The /3-model fitting results and 90% confidence 
intervals are presented in Table [31 

Assuming that the cluster is isothe rmal and in hydrostatic equilibrium, we can infer the 



mass distribution (e.g.. lNeumannll2005l ) of 



M(r) = MiL^ (4) 

where fi is taken at 0.59 for fully ionized hydrogen and helium of cosmic abundances. To 
test the validity of the isothermality assumption we fit temperatures for two portions of 
the spectral aperture: A circular region about the center of the emission and an annulus 
immediately outside the circle. The border between these sub-regions was chosen so that 
the annulus and the circle each contain roughly the same number of source counts. Results 
are consistent with the isothermal assumption but provide no additional constraint. 

We extrapolated the mass profile to obtain r^, the radius at which the spherically- 
averaged cluster density is S times the critical density at the cluster redshift, and Ms, the 
cluster mass within rs. We evaluate these parameters for 6 =200, 500 and 2500. We detect 
cluster emission above background to a radius of approximately rsoo, meaning that we see 
cluster emission from the majority of the area inside the aperture used for obtaining the 
spectral data (rsoo = 38.9" versus the spectral aperture of 49.2"). Values for Ms and rs are 
listed in Table HI We find that though the X-ray temperature, and by extension the inferred 
mass, is low co mpared to other R CS clusters of similar richness, the measured core radius 



is comparable (iHicks et al.l 120071 ). Lastly, we estimate the cluster gas mass and gas mass 
fraction at these three density contrasts by applying our /5-model for the gas density and 
comparing to our mass estimates. Results from this analysis are presented in alongside the 
mass estimates in Table HI 



Relations between observables 



3.1. Mass-temperature relation 



We compare the values of Ms calculat ed from the /j-model to the mass-temperature 
relations for M500 and M2500 determined by IVikhlinin et al.l (120061 ) in Figure [H This figure 
also includes data from six other RCS clusters, si milarly analyzed, for comparison. A full 
treatment of this larger sample will be included in lHicks et al.l (120071 ). We scaled masses by 
E{z) = H{z)/Hq to account for redshift evolution expected in the self-similar model (e.g.. 
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Bryan fc Normarull998l ). The M-Tx relation is then parametrized as 

Tx 



E{z)Ms = Mo 



5 keV 



70 ^^-^0' 



(5) 



with Mo and a taken directly from IVikhlinin et all (l2006f ): For 6 = 500, Mo = 2.93 ± 0.16 
anda = 1.61 ±0.11; for 5 = 2500, Mo = 1.28 ±0.05 and a = 1.64 ±0.06. RCS043938-2904.7 
agrees with the relation (at 90% confidence) for both 6 = 500 and 6 = 2500. 

An additional mass estimator based on X-ray data which is know n to have a tight cor- 



relation with clus ter mass is the integrated pressure Y5 = Mg^sTx (IKravtsov et al.l 12006 



Nagai et al.l 120071) . We comp a re ou r cluster data to the relation obtained through simula- 
tions done by iKravtsov et al.l (120061 ) and with the relation obtained from Chandra data by 
Nagai et al.l (120071 ). Our data agree at the 90% confidence level with both relations. 



3.2. Luminosity-temperature relation 



We adopt the luminosity-temperature parametrization of a simple power law and red- 
shift evolution following the self-similar evolution model 



E{z)-'Lx = Lo 



T 



X 



6 keV 



10^^ /ifo' erg s-\ 



(6) 



For the te mperature depende n ce exp onent and normalization, we take the local relation 
determined by lArnaud fc EvrardI (119991 ) for luminosity in the 2-10 keV band e xtrap olated 
to rsoo- This gives log^o Lq = 0.54 ± 0.03 and a = 2.88 ± 0.15. iMauehan et~aD J2006h . who 
analyzed 11 high- 2; clusters observed with XMM and/or Chandra, find this relation to be an 
acceptable fit to their data as well. Their assumed redshift evolution includes an additional 
factor to account for evolution of the density contrast needed for virialization, but they note 
that it is dominated by the self-similar evolution factor. Given the size of the errors in our 
measurements, we neglect this smaller term. We compare our extrapolated luminosity in 
the 2-10 keV band with the relation in Figure [2|, which shows that the present data are 
consistent w ith the local relatio n. Six other RCS clusters, whose complete treatment will be 
presented in iHicks et al.l (120071 ). are included in the figure for comparison. 



3.3. Mass-richness relation 



Optical richness can be quantified using the Bgc parameter introduced by lLongair &: Seldner 



(Il979l ). which has been shown to correlate well with cluster velocity dispersion, as well as 
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other observables which h ave been used for cluster mass proxies (lYee &: EUingsonll2003l ). See 



Yee fc Lopez-Cruzl (119991 ) for a detailed review of richness measures. The RCS measures 



a modified version of this parameter, -Bgcij? considering only the red- sequ ence galaxies in 
their richness calculation as an attempt to better estimate the cluster mass (iGladders &: Yee 



20051 ). The mass-richness relation is parametrized as 



Mann = lO^S^^/l + zV 10^^ M 



05 



(7) 



as m 



Gladders et all (120071 ). The measured Bg^R value for RCS043938-2904.7 of 1591±494 



(/i50^ Mpc )^-^^ (also listed i n Tab le ED was derived from the RCS data using the method de- 



scribed in lGladders fc Yed (120051 ). We use the relation obtained by lBlindert et al.l (12008. in preparation! ) 
from a dynamical study of 33 low- to moderate-redshift [z <0.6) RCS clusters. They found 
values oi A = —5.7 ± 3.4 and a = 2.1 ± 1.2 (68% confidence) and a ±0. 46 dex scatter in 
mass at a constant richness. This is consistent with the values measured by lYee fc EUingson 
( 120031 ) for the X-ray selected CNOCl sample, over a similar redshift range. The Blindert 
sample was limited in redshift and places no constraint on redshift evolution, so we take 
7 = 0. This is also cons i stent with the values and evolution inferred from the cosmological 
study of iGladders et al.l (120071 ). In Figure [3] we compare our cluster to the Blindert relation. 
Assuming that the spherically symmetric /3-model is a good representation of the system, 
RCS043938-2904.7 is a significant outlier, with far too large a richness mass for its X-ray 
mass. Optical spectroscopy measurements of RCS043938-2904.7 show evidence of significant 
substructure, as will be described in the next section. 



4. Optical spectroscopy 



Optical spectra have be en obtained for galaxies in RCS043938-2904.7 (iGilbank et al. 



20071 : iBarrientos et al.l 120041 ). The details of the observations are given in these papers. 
Briefly, the datasets comprise 34 spectroscopically confirmed members of the RCS043938- 
2904.7 system, with redshifts from LDSS-2 a nd IMACS on t he 6.5 -m Magellan Baade tele- 
scope and F0RS2 on the VLT. The spectra of iGilbank et al.l (120071 ) were classified on a scale 
of 1-5 in terms of redshift quality. For this discussion we take classes 1-3 to be secure redshifts 
and class 4 to be lower confidence redshifts. Redshift measurements were not possible for 
class 5 spectra. The VLT spectra of Barrientos et al. (2004) have been reclassified onto the 
same system. To thi s sample we add r ecent spectroscopy from F0RS2 on the VLT, similar 
to those described in IBarrientos et al.l (120041 ). These data will be presented in Barrientos et 
al., in preparation. 



As discussed in IGilbank et al.l ( 120071 ). RCS043938-2904.7 displays a clear double peak 
in the redshift histogram at 2; ~ 0.96. Figure H] shows a redshift histogram centered on 
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the mean redshift of the two components, z = 0.9558. Only galaxies within a radius of 5 
arcminutes, corresponding to a projected distance of ~ 2.0 h^Q Mpc, are considered. Fitting 
a double-gaussian to the peaks yields central locations of the two peaks which are insensitive 
to the redshift quality class under consideration. The separation of the two peaks in the 
mean rest frame ranges from 3000 km s~^ to 3250 km s~^ depending on the redshift quality 
chosen. The difference between the t wo cases is comparab le with the redshift uncertainty on 



a single measurement (~ 200 km s ^, iGilbank et al.l 120071 ). 



With the new spectroscopy, there are now enough redshifts to attempt to measure ap- 
proximate velocity dispersions for the two components. Using all galaxies of redshift quality 
4 or better gives: a = 1080 ± 320km s~^ for the z=0.9435 component and a = 560 ± 160km 
s~^ for the ^=0.9681 component, where the errors are based on jackknife uncertainties as in 



Gilbank et al.l (120071 ). Such vel ocity dispersions wou ld correspond to masses of 7.5 x 10 M 



and 1.0 X 10^^ Mq respectively (jCarlberg et al.lll997l ). However, velocity dispersion estimates 







based on so few points may be biased and the formal uncertainties are large enough that 
the values are consistent with both components having equal velocity dispersions. We note 
that both peaks in the velocity histogram contain approximately equal numbers of galaxies 
and use this, coupled with the velocity dispersion uncertainties, to adopt the conservative 
assumption that both systems have approximately equal masses. Refined mass estimates 
from dynamics will require more redshifts and such work is ongoing. Figure [5] shows the 
spatial distribution of these galaxies on the sky in a composite R, I and K band image, over- 
laid with X-ray contours. Contours begin la, assuming Poisson statistics, above the average 
number of counts per pixel in the entire image (which is essentially the background), and are 
spaced in 0.25cr intervals. Figure [6] shows the same X-ray contours along with RCS galaxy 
overdensity contours, starting at 2a and spaced in 0.25a increments, as well as the positions 
of galaxies with measured spectroscopic redshifts, including field galaxies and galaxies from 
each of the two redshift peaks. The two redshift peaks overlap in position on the sky and 
thus it is not possible to say for certain which system is responsible for the X-ray emission 
measured, or whether both systems contribute, though the latter seems most likely. From 
Figures [5] and [6] it can be seen that the peak of the X-ray emission is nearly coincident with 
the brightest cluster galaxy (BCG) of the z=0.94 cluster, and that the emission is elongated 
in the direction of the 2=0.96 cluster's BCG. This adds further support to the idea that 
the observed X-ray emission is likely due to a combination of emission from each of the two 
systems. 

Figure [31 as noted above, shows Bg^R plotted against M200 for our RCS043938-2904.7 
and a selection of other RCS clusters. RCS043938-2904.7 lies well outside the region bounded 
by the 1-a scatter, appearing significantly richer than its X-ray mass would suggest. This is 
not surprising given that the spectroscopy discussed above suggests that RCS043938-2904.7 
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is an extended structure containing two comparable, less rich systems. Since these systems 
lie so close together in redshift that their red sequences overlap, the red-sequence Bg^R 
measurement is the sum of that of both these systems. A more reasonable estimate of the 
true richness associated with the measured mass X-ray mass, in light of the double-peaked 
velocity distribution, would be to halve the measured BgcR. Doing so brings the X-ray mass 
and the richness mass into agreement, within the large scatter. 



5. Discussion 

RCS043938-2904.7 has an inferred X-ray mass which is one to two orders of magnitude 
lower than the mass indicated by its red-sequence richness, assuming that the X-ray emission 
is from a single, spherical, isothermal gas distribution and that the richness is associated with 
a single cluster (Msoo.x = 4.6 i?;? xlO^^ versus M2oo,B,,n = 1-1 ^ 10^^ ± 0-46 dex M©). 
This is a significant discrepancy for two reasons. First, as mentioned before, richness is 
well correlated on average with other observable properties. Secondly, clusters with masses 
above 10^^ Mq are quite rare, especially at 2; ~ 1. This means that if the richness mass 
estimate of RCS043938-2904.7 were truly indicative of its size, it would be among the largest 
overdensities in the observable universe, while also being very underluminous. Yet the X-ray 
properties, which suggest a much less massive object, are consistent with the expected M5- 
Tx and Lx-Tx relations, meaning that the temperature of the plasma is consistent with its 
observed distribution. One possible interpretation of these results is that the red-sequence 
galaxies used to measure the optical richness and the X-ray emitting gas trace different 
volumes. Specifically, while the X-ray luminous plasma is confined within one or more deep 
gravitational potential wells, the galaxy population may occupy a more extended region 
which is not yet dynamically relaxed. The existence of such a structure in RCS043938- 
2904.7 is supported by its two-peaked radial velocity distribution. 

Merging and dynamically ac tive clusters are exp ected and observed to be increasingly 



common at higher redshifts (e.g.. iJeltema et al.ll2005l ). Thus any cluster sample will include 



an increasing fraction of dynamically young systems at higher redshift. Furthermore, N-body 
simulations indicate that a significant fraction of cluster samples selected using broadband 
color discrimination will be systems whose galaxy members are not predominantly associated 
with a single, large potential well, but rather are spread amongst a number of smaller, but sti ll 



physically associated, dark matter halos distributed along the line-of-sight (jCohn et al.ll2007l ). 
Simulations of the RCS technique on mock catalogs tuned to reproduce observables such as 
the observed galaxy color distribution and the two-point correlation function sho w that we 



expect false-positive cluster detections to occur at a frequency of ~5% in the RCS (iGladders 



2002 ) . This agrees with i nitial results based on a small number of clusters (iGilbank et al, 
2OO7I : Isiindert et alj[2007h . 



Since the red-sequence galaxies used in the RCS richness measurements form very early 
{z > 2) in high density regions, overdensities of red-sequence galaxies can be associated 
with very large structures which may not be relaxed during the observed epoch. The X-ray 
luminous plasma, on the other hand, will be confined to one or more gravitational potential 
wells whose size is limited by the collapse timescale. In the hierarchical collapse paradigm of 
a ACDM universe, the large, virialized wells seen in local clusters do not develop until long 
after the galaxies have formed. In this scenario, for any evolving cluster, there is likely to 
be an epoch at which the apparent richness overestimates the virialized mass. 

Moreover, the misinterpretation of the X-ray emission from a complex, dynamically 
young cluster as a single, virialized structure can lead to an overestimate of the inferred gas 
mass fraction and an underestimate of the total mass, which can be seen as follows. Let us 
consider two models for the X-ray emission seen in RCS043938-2904.7. Our one-component 
model (Model I) is the set of assumptions used in our X-ray analysis thus far: X-ray emission 
is from a single, spherically-symmetric distribution of gas in hydrostatic equilibrium with 
a gravitational potential which is well- described by a /5-model. The two-component model 
(Model II) assumes that there are two nominally identical (i.e. with the same values of (3, tq 
and Tx), virialized gas distributions, separated along the line of sight, and that each of the 
two gas distributions is responsible for half of the measured X-ray flux. Each component is 
assumed to obey the local M^-Tx and Lx-Tx relation. 

Because it assumes the flux is halved between the two components, the luminosity 
inferred per cluster in Model II is half that of Model I. The luminosity scales as the square 
of the gas mass, so since the spatial distribution is the same for each cluster, the inferred 
gas mass and thus also the electron density will be reduced by a factor of a/2 in each of the 
clusters relative to Model I. Since our total mass estimate for each cluster depends only on 
the temperature of the X-ray emitting gas and its distribution in the plane of the sky, each 
of the two components in Model II will have the same total mass as the single component in 
Model I. This means that the inferred gas mass fraction will be a factor of \/2 lower using 
Model II versus Model I. The total mass, total gas mass and gas mass fraction results using 
each of the two models can be found in Table |H 

Similarly, n identical components along the line of sight would reduce the inferred gas 
mass per component and the overall gas mass fraction by a factor of ^/n in addition to reduc- 
ing the luminosity per component by a factor of n. In principle, luminosity measurements 
and the Lx-Tx relation constrain the number of components allowed, but our luminosity 
errors are too large to distinguish between Model I and Model II. 
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Mathiesen fc EvrardI (120011 ) . using a suite of hydro dynamic cluster simulations, found 
that during merger events the ICM spectral fit temperature will underestimate the mass- 
weighted ICM temperature by ~ 20%, because the cool, denser inflowing gas will dominate 
the emission over the gas already heated by the merger. This means that for a dynamically 
young system, such as RCS043938-2904.7, where even the most virialized components are 
likely to have undergone recent mergers. X-ray mass estimates for those clusters may still 
be below the virialized masses by sever al tens of percent. In a recent detailed study of 
the cluster CI 0024-1-17, |Jee et al.l (120071 ) found a similar underestimation of the true mass 
distribution due to the assumption of a single virialized mass structure rather than two 
components extended along the line of sight and in a state of ongoing dynamic interaction. 

In summary, the assumption of spherical symmetry and virial equilibrium for a cluster 
system containing extended line-of-sight structure and dynamic evolution can lead to over- 
estimation of the gas mass fraction and underestimation of the total mass. The baryon mass 



fract i on in relaxed galaxy c lusters is expected to be a universal quantity (e.g., IWhite et al. 



19931 : IVikhlinin et al.ll2003l ). suggesting that the comparison of cluster gas mass fractions 
inferred by assuming spherical symmetry to canonica l cluster valu e s mig ht be used to infer 
the presence of line-of-sight structure. As noted by IWhite et al.l (119931 ) . the cosmological 
baryon mass ratio, r2b/f2m5 provides an upper limit to the total gas r aass fraction. This r atio 
from the WMAP three-year data results is Vt^/VL^ = 0.175 ± 0.012 JSpergel et al.ll2007h . 



Galaxy cluster samples at high redshift are expected to have lower cluster ga s mass 
fractions than local samples (IHicks et al.l 120071 : ISadat et al.ll2005l : lEttori et al.ll2004l . e.g.,). 
Also, a positive correlation between temperature and gas mass fraction has been observed 



(e.g. . IVikhlinin et al.ll2006l : ISanderson et al.ll2003l ). This suggests that a high-redshift cluster 
with a low X-ray temperature, such as RCS043938-2904.7, ought to have a correspondingly 
low gas mass fraction. Instead, we observe a very high gas mass fraction assuming that it is 
a single, spherical matter distribution, especially measured within a large radius. At t^qq, we 
find /gas, 500 = 0.17±q;o8. We observe cluster emission to approximately this radius, so this gas 
mass fraction estimate is unlikely to have the extrapolation errors that the measurement at 
r2oo might. Assuming that Model II is more appropriate for the physical state of RCS043938- 
2904.7 lowers the gas mass fraction towards what is expected. 

Would this pair of clusters be physically associated? As an order-of-magnitude argu- 
ment, we note that the richness mass (~lxlO^^ M©) is roughly equal to that of a sphere 
of radius ~13 Mpc with the cosmic density of z = 0.96. Optical spectroscopy shows that 
the redshift separation of the two components in RCS043938-2904.7 is Az ~ 0.005, or a 
physical separation of about 12 Mpc in the Hubble flow. This rough agreement supports 
the conclusion that RCS043938-2904.7 is an incompletely-virialized system which is still ap- 



-13- 



proaching equilibrium, and that the galaxy distribution traces unvirialized matter extended 
along the line-of-sight, as well as the virialized matter traced by the X-ray gas. The free-fall 
collapse time for the two components given by the richness mass is approximately equal 
to the lookback time, indicating that this cluster would be nearly virialized by about the 
present epoch. 

Additional observations are required to determine the physical state of RCS043938- 
2904.7. Due to the low number of source photons from this cluster, there are large errors in 
the X-ray luminosity and temperature. Improvements on both would increase the precision 
of the X-ray mass estimates and the gas mass fractions. An independent mass estimate, 
such as might be obtained from weak Icnsing, could also help distinguish which mass esti- 
mate (richness or X-ray) is more appropriate, as well as giving additional insight into the 
distribution of matter. A weak lensing measurement, though difficult, would be particularly 
interesting for this system if it could trace the spatial distribution of cluster mass perpen- 
dicular to the line of sight, possibly revealing more substructure. Additional spectroscopy, 
currently underway, will help to further constrain line-of-sight substructure. 



6. Summary 

We present new CXO observations and optical spectroscopy measurements for the high- 
redshift galaxy cluster RCS043938-2904.7 detected in the Red-sequence Cluster Survey. We 
have obtained a spectral temperature and estimated total mass, gas mass and gas mass 
fractions from the X-ray data. The measured X-ray properties and mass estimates agree 
well, within errors, to the locally derived Lx-Tx and M^-Tx relations for galaxy clusters, 
but not with the local M2oo-BycR relation. Spectroscopy of RCS043938-2904.7 shows a 
two-peaked velocity distribution, suggesting substructure along the line-of-sight. Assuming 
that there is a single, viriahzed mass distribution in RCS043938-2904.7 results in gas mass 
fraction measurements which are higher than others in the RCS, which may also indicate 
line-of-sight structure. Accounting for the substructure with a two-component model for the 
mass distribution brings the richness mass into agreement with the X-ray mass and the gas 
mass fractions into agreement with expected values. Weak lensing measurements, additional 
X-ray observations and additional spectroscopy would be useful in further constraining both 
the total amount of mass in the system and its distribution. 
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8. Tables and Fig 


^ures 


Obsids 


Observation Dates" 


Livetime^ Counts'^ 


3577, 4438 


04/16/2003, 06/06/2003 


93.3/105.0 359/2093 



Table 1: Observation Log. a: Dates are in m/d/y format, b: Livetime is in ks, and 
split into good/total amounts, c: Counts are in the 0.3-7.0 keV band, and split into 
source/source+background totals. 



Redsliift 


rp a 


T b 


r b 


^bol 




0.9558 


1 5+1-0 


4.0±|i 


2.l±\l 


6.9±|:° 


13.7±?i 



Table 2: Spectral Analysis Results, a: Temperature is in keV. b: Luminosities in 10^^ erg 
s~i. Subscripts refer to the rest frame energy band; the bolometric luminosity is estimated 
by the 0.05-50 keV band, c: Measured flux in the 0.3-7.0 keV band, in units of 10"^^ erg 
cm~^ s~i. All values are from within the 49.2" radius spectral aperture. 



-'0 




ro " 


13 BgcR 


3 7+2-5 

"J- ' ="=2.0 


3.9 + 1.5 


14.7+1.9/117+15 


0.56+°:°^ 1591+464 



Table 3: /?-model Results and Optical Richness, a: Central surface brightness, in 10" erg 
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cm ' s arcmin b: Central electron number density, in 10 cm c: /3-model core 
radius, in arcseconds/ kpc. d: Optical richness, in {K^q Mpc)^ ''^, as measured by the RCS 
jciadders fc Yeelboosh . 
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Model I 






Model II 




6 


rs " 


Ms' 




f 

J gas,S 






f 

J gas,S 


200 


510±^^° 


4.6±f? 


1.05±°jl 


0.23±H? 


9.2±i2^o 






500 






0.45±°:?^ 


n 1 7_(_0.07 


5 2+'^-'^ 

^■^=■=2.0 


0.64±0;i 




2500 


90±tl 


U.OOICq 3Q 




0Q+°-°^ 


"J- ' ="=0.6 


0.04±0;J1 


Q 07-^0.02 

="=0.03 



Table 4: Total Mass, Total Gas Mass and Gas Mass Fraction Estimates. Mass values are 
for a single component in Model I and for the sum of two components in Model II. a: The 
physical radii corresponding to each of the three density contrast radii, in hjQ kpc. b: Total 
mass for each model and density contrast, in 10^^ h^Q Mq. c: Gas mass for each model 
and density contrast, in 10^^ ^m^'^ ^o- d: Gas mass fractions for each model and density 
contrast, in units of h^Q^^- 




Fig. 1. — Mass- Temperature Rela tion for M25nn ( a t left ) and M500 (at right). The solid 



lines are the relation obtained by IVikhlinin et al.l (120061 ). The red square point on each 



plot indicates RCS043938-2904.7. The blue triangle points indicate a selection of other RCS 
clusters which were similarly analyzed and are provided only for comparison. A detailed 
analysis of the set of RCS cl usters and their red shift evolution is beyond the scope of this 
paper and will be provided in lHicks et al.l (120071 ). 



- 19 - 



1 - 



0.1 - 




Tx (keV) 



Fig. 2. — Luminosity- Temperature Relation. Luminosity is in the 2-10 keV band, extrapo- 
lated to the radius r2oo- Redshift evolution is scaled out according to the self-similar model. 
The red square point indicates RCS043938-2904.7. The blue triangle points indicate a selec- 
tion of other RCS clusters which were similarly analyzed and are provided only for compari- 
son. A detailed analysis of the set of RCS c lusters and their re dshift evolution is beyond the 
scope of this paper and will be provided in iHicks et al.l ( 2007 ). 
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gcR 



1000 



1500 



2000 



Mpc) 



1.77 



Fig. 3. — The M2oo^BgcR relation, plotting X-ray mass against richness. The solid line is 
the Blindert et al. relation and the dotted lines above and below are the la scatter (±0.46 
dex). The red square point indicates RCS043938-2904.7 assuming that all the richness is 
associated with a single component (Model 1). The green diamond point indicates each of 
the two equal- mass components within RCS043938-2904.7 assuming Model 11, with half 
the observed richness associated with each component. The blue triangle points indicate 
a selection of other RCS clusters which were similarly analyzed and are provided only for 
comparison. A detailed analysis of the set of RCS cl usters and their re dshift evolution is 
beyond the scope of this paper and will be provided in lHicks et al.l (120071 ). 
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Fig. 4. — Rest-frame velocity histogram centered around the mean velocity of the 
RCS043938-2904.7, z=0.9558. Hatched histogram shows secure redshifts and open histogram 
shows lower confidence measurements. The overplotted curves show best double-gaussian 
fits to these, solid and dashed curves respectively. Uncertainties arc sufficiently large that 
velocity dispersion measurements for each of the two components are consistent with equal 
velocity dispersions. 
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Fig. 5. — Three-color image (Red, green and blue are K, I and R bands, respectively) of 
RCS043938-2904.7 with overlaid green X-ray contours. North is up and east is to the left, and 
the image is 2' (797 kpc) on a side. The X-ray image was smoothed using a 0.75" gaussian 
kernel and X-ray point sources were removed using the CIAO script dmf ilth. Contours 
begin Icr above the mean counts per pixel (assuming Poisson statistics) and contours are 
separated by a 0.25cr spacing. 
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Fig. 6. — In red, X-ray contours as in Figure [5] and in blue RCS galaxy overdensity contours, 
starting at 2cr and increasing in 0.25a increments. Points indicate galaxies with measured 
redshifts from optical spectroscopy. X points indicate galaxies with which are not within a 
±5000 km s~^ window of the mean cluster redshift, box points indicate galaxies associated 
with the 2=0.94 cluster, and diamond points indicate galaxies associated with the 2=0.96 
cluster. 



